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Ph I Abstract: We consider a variant of the classic Steklov eigenvalue problem, which 

■ arises in the study of the best trace constant for functions in Sobolev space. We 

^ ■ prove that the elementary symmetric functions of the eigenvalues depend real- 

d '• analytically upon variation of the underlying domain and we compute the cor- 

responding Hadamard-type formulas for the shape derivatives. We also consider 
isovolumetric and isoperimetric domain perturbations and we characterize the 
corresponding critical domains in terms of appropriate overdetermined systems. 
Finally, we prove that balls are critical domains for the elementary symmetric 



. functions of the eigenvalues subject to volume or perimeter constraint. 
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1 Introduction 

Let N >2 and f2 be a bounded domain {i.e., a bounded connected open set) in 
of class C^. We consider the eigenvalue problem 

Au = u, in Q, 

(1-1) 

g = \u, on dn, 
and we represent its eigenvalues by means of a divergent sequence 

< Xi[n] < X2[n] <■■■< Xj[n] < ... 

where each eigenvalue is repeated according to multiplicity. Here u denotes the 
outer unit normal to dQ. In this paper we study the dependence of Xj[fl] on Q. 

After the pioneering paper by Steklov [16], the homogeneous boundary con- 
dition du/dv = Xu is called Steklov boundary condition and it appears in the 
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study of many boundary value problems. In the literature, the Steklov condition 
is often imposed to harmonic functions in Q and the eigenvalues of the problem 



Au = 0, in f2, 
p = Am, on dn, 



(1.2) 



are usually called Steklov eigenvalues. Problem (11. 2p has important physical 
meaning and many properties of the Steklov eigenvalues are known. In particular, 
for N = 2 problem (11.21) describes the vibration of a free membrane the total mass 
of which is uniformly distributed along the boundary. On the other hand, problem 
(11.11) seems to find its natural motivation in the frame of the Theory of Sobolev 
Spaces and in particular in connection with the problem of traces. Recall that 
functions of the Sobolev space W^''^{fl) admit a trace in L'^{dfl) and there exists 
a constant C such that 



for all u E W^'^(ri). Moreover, the trace operator from W^''^{il) to LF'{dVL) is 
compact. By simple considerations, it turns out that 



\i\VL\ = min 

Trn^O 





\Vu\ 




\u\ 




Idfl 


\u\ 





and that the minimizers are exactly the eigenfunctions corresponding to Ai[17]. 
In fact, problem (II. ip is a formulation of the Euler-Lagrange equation associated 
with (II. 4p . Note that Xi[Q] is the reciprocal of the best Sobolev trace constant 
C for which (I1.3P holds. 

In the last decade problem (II. ip and its nonlinear generalizations have at- 
tracted the attention of several authors, see, e.g., the rather widely quoted paper 
by Martinez and Rossi [H]. In particular, Rossi flS] has conjectured that the the 
ball maximizes Xi[Q] among all domains with fixed volume (this is known to be 
true for the first nontrivial eigenvalue of (II. 2p . see Weinstock [17] and Brock [2], 
see also Henrot [7]). In order to support his conjecture, Rossi [15] has computed 
the shape derivative of Xi[^l] and has proved that the ball is a critical point for 
Ai[0] under volume constraint. 

In this paper, we develop the results of [15] and we consider all eigenvalues 
Aj[n]. The main difficulty in dealing with higher eigenvalues is related to their 
multiplicity which may change when the domain is perturbed. This leads to 
complicated bifurcation phenomena, which clearly do not occur in the case of the 
simple eigenvalue Ai[f2]. The investigations carried out in the papers [HI [121 E], 
have pointed out that in the case of multiple eigenvalues it is natural to consider 
the elementary symmetric functions of the eigenvalues, which typically depend 
real-analytically on the parameters involved in the problem. In this paper, we 
adopt the approach of [12]. Namely, we consider families of domains 4>{fl) param- 
eterized by suitable diffeomorphisms (p of class defined on the fixed reference 
domain Q and we prove that the elementary symmetric functions of the eigen- 
values of problem (II. ip on (f){Q) depend real-analytically on 0, see Theorem 13.21 
Moreover, in formula (13. 4p we compute the Frechet derivatives of such functions 




(1.3) 
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and we characterize the critical transformations (f) subject to the volume con- 
straint Vol(i7) = const or the perimeter constraint Per(0(f2)) = const, see Theo- 
rem 14.81 This leads to the formulation of the appropriate overdetermined systems 
(I4.13p . fl4.14p . It turns out that if 4>{fl) is a ball, then is a critical point for 
the elementary symmetric functions of the eigenvalues with volume or perimeter 
constraint, and the corresponding overdetermined systems are satisfied. 

We note that in the terminology of domain perturbation theory (c/. Henry [6]), 
in this paper we adopt the 'Lagrangian' point of view, whilst the calculations in 
[15] are performed in the 'Eulerian' form. 

2 Notation and preliminaries 

In this paper the elements of are thought as row vectors. The inverse of a 
matrix A is denoted by as opposed to the the inverse of a function / which is 
denoted by f^~^\ The transpose of a matrix A is denoted by A"^ , and the inverse 
of the transpose of A is denoted by A'"^. The product of two matrices A and B 
is denoted by AB. According to this notation, W1W2 is the scalar product of two 
vectors wi,W2 G M^. The symbol V denotes either the gradient of a real-valued 
function or the Jacobian matrix of a vector- valued function. The differential of a 
function F at a point a is denoted by dF(a) and the value of dF(a) computed in 
b is denoted by dF{a)[b]. 

Let Q he a. bounded domain in of class C^. Here and in the sequel it is 
understood that N >2. We consider the weak formulation of problem (II. ip 



in the unknowns u G W^''^{Q) (the eigenfunction) , A G M (the eigenvalue). Here 
W^'^{n) denotes the Sobolev space of those functions in L^(f2) with first order 
weak derivatives in L^(f2) endowed with its standard norm (functions in L^(fi) 
are assumed to be real- valued) and da denotes the (A^ — l)-dimensional surface 
measure in dfl. Note that the integral in the right-hand side of fl2.ip is well- 
defined since functions in the Sobolev space admit a trace in L'^{dfl) in the sense 
that there exists a continuous linear operator Tr of W^''^{fl) to L'^{dfl) such that 
Tiu = u\QQ for all u G W^''^{VL) fl C{VL). Occasionally in the sequel we shall also 
use other standard Sobolev spaces iy™'^(f2) of functions with weak derivatives 
up to order m in U'{VL)] we refer to Burenkov [5] for the basic properties of these 
spaces. 

We now recall a standard procedure that enables us to reduce the study of 
problem (12. ip to the study of an eigenvalue problem for a compact self-adjoint 
operator in the Hilbert space W^^'^iyL) (which is equipped with its standard scalar 
product defined by the left-hand side of (12. ip ). 

We consider the Laplace operator A as an operator from W^'^{VL) to its dual 
(iyi'2(n))' defined by A^M = - j^VuV^'^dx for all u, G W^^'^{9) and we 
recall that the operator / — A is an isometry of W^''^{VL) onto {W^'''^{VL))' . Here 
/ is the natural embedding defined by /[tt][v2] = J^u(pdx for all u,{p E W^'^{Q). 
Moreover, we consider the embedding J of L'^{dQ) into {W^'^{fl))' defined by 
J[m][v9] = Jg^ULpda for all u G L^{dQ), ip G W^'^{Q). Clearly, equation f l2.ip is 




(2.1) 
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equivalent to equation 



{I - A)-^ oJoTt[u] = X-^u (2.2) 

in the unknown u G W^''^{Q). Note that by using (p = u as a test function in 
(12.11) . it is immediate to see that any eigenvalue A is positive. Thus, we have the 
following 

Lemma 2.3 The operator Sn = {I — A)~^ o J o Tr is a compact non-negative 
self-adjoint operator in the Hilbert space W^''^{Q). Moreover, S^u = fiu for 
u G W^''^{VL), /i G M \ {0} if and only if u is an eigenfunction of problem \2. 1\) 
corresponding to the eigenvalue A = 

Proof. By equality 

< SnUi,U2 >H/i.2(n)= {I - A)[(/ - A)"^ o J o Tr [mi]][m2] = J o Tr [mi][m2], 

for all ui,U2 G W^''^{Q) and by the symmetry of the pairing J o Tr [■] [■] we imme- 
diately deduce that Sq is a non-negative self-adjoint operator. The compactness 
of Sq follows by the compactness of the trace operator Tr (see, e.g., Kufner, John, 
Fuci'k jHl Thm. 6.10.5]). The proof of the last part of the statement is straight- 
forward. □ 

It is easy to see that KerSn = Wq''^{Q), where Wq''^{Q) is the closure in 
W^''^{Q) of the space of functions of class C°° with compact support in Q. Since 
the operator Sq is compact and self-adjoint in the Hilbert space W^''^{Q) and its 
kernel has infinite codimension in W^''^{Q), the spectrum of ct(S'q) is discrete and 
ciSn) \ {0} consists of eigenvalues of finite multiplicity which can be represented 
as a non-increasing sequence /ij[S'n], j G N converging to zero (as customary, 
we repeat each eigenvalue as many times as its multiplicity). Thus we have the 
following 

Corollary 2.4 Let Q be a bounded domain in of class . The eigenvalues 
of problem i\2. 1\) are positive, have finite multiplicity and can be represented as a 
non- decreasing divergent sequence Xj[Q], j G N where each eigenvalue is repeated 
according to its multiplicity. Moreover, 



XjlQ] = min max -^^ — j; — — — , (2.5) 

IT 7/3^(1 Jail I I 

EnWl'^(Q,)={Q} "^'^ 

dim_B=j 

for all j G N and formula \1.4\j holds. Furthermore, Ai[i7] is simple and its 
eigenfunctions do not change sign in Vt. 

Proof. The proof of the first part of the statement can be easily deduced by 
Lemma 123] and standard spectral theory and by recalling that Aj[fi] = /ij^[5'n] 
for all j G N. For the second part we proceed exactly as in the well-known case 
of the Dirichlet Laplacian. Namely, let m be a nonzero eigenfunction associated 
with Ai[f2]. Since m is a minimizer in (II. 4p then also \u\ is a minimizer in (II. 4p . 
Thus, \u\ is an eigenfunction associate with \i\Vt\ and by the Harnack inequal- 
ity \u\ cannot vanish in VL. Hence u does not change sign in VL. Finally, given 
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two nonzero eigenfunctions Ui,U2 associated with Xi[Q] and c G M such that 
J^ui — cu2dx = 0, it follows that the eigenfunction ui — CU2 must be identically 
zero, hence ui = CU2- □ 



Remark 2.6 Another way of reducing problem ^2. ij) to an eigenvalue problem 
for a compact non-negative self-adjoint operator in Hilbert space is to consider 
the operator Tr o (J — A)~^ o J in L'^{dQ). This operator may be used instead 
of Sfi in all our arguments below without any essential changes. With regard 
to this, we note that, in order to normalize eigenfunctions, the scalar product 
of L'^{dQ) would be more natural than the scalar product of W^'^{Q) (see also 
Remark \3.5\ in Sectionl^. However, the relation between the eigenfunctions of 
Tr o (/ — A)"-"^ o J and those of problem Ii2.1\) is slightly involved. Thus, we prefer 
to use the operator Sq since its eigenfunctions are functions defined on fl and 
coincide with the eigenfunctions of problem Ii2.1\) . 

3 Domain perturbation 

Let f2 be a fixed bounded domain in of class C^. We consider problem (12. ip on 
a class of domains diffeomorphic to Q. Namely, we consider the class of domains 
{(p{fl) : (p G An} where An is defined by 

An = {(l)^C^i^; K^) : (j) is injective and detV0(x) ^ 0, for all x E Cl} , 

(3.1) 

and C^(r2;]R^) is the space of the functions from Q to of class C^. Note 
that if G An then (f){Q) is a bounded domain of class in and that 
d(p{fl) = 0(c}f2). Thus, it is possible to set the eigenvalue problem (12. ip on 0(fi) 
for all G An and to consider the eigenvalues Aj[0(i7)] as functions of 0. For 
simplicity, we set Aj[0] = Xj[(j){fl)]. 

In this section, we study the dependence of upon variation of G An- It 
is understood that the space C^(f2 ; M^) is endowed with the C^-norm defined by 
U\\ = Eti E\a\<2 P"0dloo for all = (01, . . . , 0^) G C\n ; M^), where || ■ |U 
is the standard sup-norm. It turns out that the set An is an open set in the space 
C^(r2;M^) (see jTOl Lemma 2.2]). Thus, it makes sense to study differentiability 
and analyticity properties of the maps Xj[-] defined on An- 

Following [TT], we fix a finite set of indexes F G N and we consider those 
maps G An for which the eigenvalues with indexes in F do not coincide with 
eigenvalues with indexes not in F; namely we set 

An[F] = {<peAn-- A,[0] A;[0], V j G F, / G N \ F} 

and 

BnlF] = {0 G An[F] : A,J0] = A,,[0], V j,,j2 G F} . 

Then we can prove the following theorem, where Vrviy) denotes the orthog- 
onal projection of Vv{y) to the tangent plane of 90(f2) at the point y G d(f){Q). 
Note that by standard elliptic theory (c/. Agmon, Doughs and Nirenberg p]; see 
also the recent paper by Castro Triana [1]), any eigenfunction v in 0(fi) is of class 
Vr^'^(0(f2)) hence Vv admits a trace in L^(90(f2)). Thus, Vr^ and the normal 
derivative dv/dv are well-defined and |Vf p = iV^f P + \dv/dv\^. Although it is 
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not necessary, it may be useful to recall that C^'°-type results (c/., e.g., Lieber- 
man [13]) imply that the gradient of v is in fact continuous up to the boundary 
of 4>{^), hence Vtv and dv/dv are actually defined everywhere in 

Theorem 3.2 Let Q be a bounded domain in of class and F a finite subset 

ofN. Then An[F] is an open set in (7^(0 and the elementary symmetric 
functions 

Af,/.[0]= Y1 h=l,...,\F\ (3.3) 

ji<-<jh 

are real-analytic in Aq\F]. Moreover, if (p & Qci[F] is such that the eigenvalues 
Xj[4>] assume the common value Xpy^] for all j G F, then the differential of the 
functions Ap^h the point is given by the formula 

+ (1 - XF[(t>]H - Xl[(t>])vf] o (t>(-^^)u^da, (3.4) 

for all ip G (7^(0 ;M^), where {vi}i(zf is an orthonormal basis in W^''^{(I){Q)) of 
the eigenspace corresponding to Xf[4>], and H = divu is the mean curvature of 
90(0) (the sum of the principal curvatures). 

Remark 3.5 By looking at the weak formulation h2. it is immediate to real- 
ize that if one considers eig en functions vi normalized with respect to the scalar 
product of L'^{d(f){fl)) then the factor Xp[(f)] in formula Iji3.4^ should be replaced 
by Xp~^[(f)]. In particular it follows that formula ^3.4^ computed for \F\ = h = 1 
agrees (up to a misprint) with the formula in Rossi (TSj (1.5)] concerning the first 
eigenvalue. 

The rest of this Section is devoted to proving Theorem 13.21 We note that the 
proof of Theorem 13.21 is based on a general result in [11] concerning families of 
compact self-adjoint operators in Hilbert space with variable scalar product. In 
order to apply that result we need to consider the eigenvalue problem fl2.ip on 
the variable domain (j){VL) and pull it back to the fixed domain VL by means of a 
change of variables. In particular, for a fixed cf) G An-, it is convenient to consider 
the operators A, J, and J on (j){VL) defined in Section 2, and pull them back to 
Vt. In this way we obtain the operators A^, of W^''^{Vt) to (W^''^{VL))' and 
of L2(aO) to (l^i'2(0))' defined as follows: 

A^ [u] M = A [^i o 0(-i)] [p o 0(-i)] , I^[uM = l[uo 0(-i)] [p o 0(-i)] , (3.6) 

for all u,ip e ly^'^(O), and 

J^MM = J[no0(-i)] [^o0(-i)], (3.7) 

for all u G L^((90), G W^''^{VL). By means of standard calculus and a change of 
variables one can check that 

A,^ [u] [ip] = - f Vu (V0)"^ (V0)"^ V<^^ |detV0| dx, (3.8) 
Jo. 
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I^[u][(p]= / uip\deW(f)\dx, (3.9) 
Jn 

for all u,!f e W^'^{Q), and 

J^[u][ip]= / |z/(V0)"^| |detV0|rfa, (3.10) 

Formulas (13. 8p and (13.91) suggest introducing in W^''^{Q) the following scalar 
product 

<Mi,M2>0= / (Vmi(V0)~^(V0)"^ Vm^ + MiMs) |detV0|c?x, (3.11) 
Jn 

for all Mi,M2 e ly^'^(fi). We denote by W^'^{Q) the space iy^'^(f2) equipped 
with the scalar product defined by (13. lip . Note that the scalar product (13. lip is 
equivalent to the standard scalar product of W'^''^{Q), hence VF^'^(fi) is in fact a 
Hilbert space. Note also that the map Q of W^'^{(j){Q)) to W^'^i^Q) defined by 
C^[v] = V o (f) for all v G W^''^ is an isometry. By setting 

= (/^ - A^y^ o J^o Tr, 

it is immediate to verify the validity of the following 

Lemma 3.12 The operator T^p ofW^''^{^l) to itself satisfies 

hence it is unitarily equivalent to S^{yi). In particular, is a compact non- 
negative self-adjoint operator in W^''^{Vl) the eigenvalues of which coincide with 
the eigenvalues of S^q). Moreover, v G W^''^{(f){fl)) is an eigenfunction of S^q) 
if and only if v o (p is an eigenfunction ofT^. 

In order to prove formula (13. 4p we need the following technical lemma. To 
avoid heavy notation, we use the same symbol u to denote the outer unit normal 
to the boundaries of different domains (the context is usually clear enough to 
understand which domain u refers to; for example, u refers to dQ in (I3.14p and 
d(j){n) in IKm ). 

Lemma 3.13 Let Q be a bounded domain of class . Let u G W'^'^{VL) be fixed. 
The function B of An to M defined by 

B[(t)\= f M |z/(V0)~^| |detV0|rfa, (3.14) 
Jan 

for all (j) G An, is real- analytic and the differential at the point cf) G Aq, is given 
by 

dB{(P)[i)]= f (Hv + ^-\{i)0(P^-^^)v'^da- f (7/^o0(-i))Vt;^rfa, (3.15) 
for allipe C\n;R^), where v = mo0(-i). 
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Proof Since detV<^ does not vanish on dQ, the weight |z/(V</))~^| |detV<^| 
does not vanish on dQ. Thus, by a continuity argument infgn I'^lV^)^^] |detV0| 
is strictly positive in a suitable neighborhood in An of any fixed G An- Thus, 
the map B is real-analytic since it is the composition of real-analytic maps. 

By standard calculus, it is easy to see that differentiating detV0 with respect 
to yields 

[(d (detV0) o detV0("^) = div {ip o , (3.16) 

for all ipeC^in; M^). Moreover 

Hi,.v7^wii_ KV0)-Mv^(V0)-^ + (V0)-^(VV^)^](V0)-^^^ 

'^''^^^'^^ 2|z/(V0)-i| ■ ^ ^ 

We note that the outer unit normal to d(f){Q) at the point 0(x) G d(j){Q) is given 
up to the sign by 

Kx)(V0(3;))-^ 

|z/(a;)(V0(a;))-ir ^ ' ^ 

where z/(x) is the outer unit normal to dQ at the point x G Thus by combining 
fl3.17p . f l3.18p and by means of a change of variable (see also formula f l3.10p ) and 
by the Divergence Theorem applied in 0(f2) we have that 



ud\u{V(j))-^\['^]\deW(j)\da 



an 



I- [ t;z/[V(^o0(-i)) + V(^o0(-^))^]z/^da 
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= - [ div(t;V(V'o0(-iV^)^2/- (3-19) 

J<j>{Q) 

In fl3.19p and in the sequel it is understood that the unit vector field z/(?/) defined 
for all y G d(j){fl) is extended as a unit vector field of class in a neighborhood 
of d(j){Q) and then extended in any manner as a vector field of class defined 
in the whole of 0(f2) (see also Henry [BJ p. 15]). By fl3.16p and fl3.19p it follows 
that 

dB[^][ij]= [ vdiY{ij o (f)^'^^)d(T - [ div(t;V(V^ o 0(-i))z/^)rf?/ 



I t;div(^ o 0(-i))rfa - / Vf V(^ o 0(-i))z/^rf?/ 
Jd^{n) J(i>(n) 

- [ t;div(V(V^o0(-i))z/^)c/?/. (3.20) 
J(h(n) 



To proceed with these computations we need to differentiate the vector field u 
twice. Thus, since u is of class C^, we approximate z/ by a sequence of vector 
fields z/(") = (z/J"\ . . . , i/Jj"^), n G N such that ^ z/ in Ci(0(fi)) as n ^ oo. To 
shorten our notation we set ( = ipo(f)^~^\ By applying repeatedly the Divergence 
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Theorem, we obtain 

^iJ^m oyrdys dyr dys Jd^(n) 

^^iJnn) oyrdys j^iJdH^) 9ys 

+ E/ + (3.21) 

r^lJ<f>m "Vs oyrdys 

and 

^(n) rl^J^i^) ^y^^yr 9ys dyr 



Thus 



f ( in) T.r ^ f 9V dCr („) , dCrdu^ 

+ vCri^s^^da - Y] / ' +Cr o dy 

r^iJd<p{n) oyr -^Z/r^l/s dyr dys 

Jd<i>{n) J(j){n) oys oyr oyr oys 

r^iJdH^) oyr rT^i-^ 't'i^) ^yr^Vs dyrdys 

-TT^s'-^^dy- [ {u^-Wv^Xu^da. (3.22) 
dVr dys Ja^(n) 



I Vv\/C{u^"'Ydy+ I i;div(VC(z/^"^)^)rfy = / (z/(")z/^)t;divCrfa 

J<i>(Q.) Jd(j){n) 

Jd^{Q) Jd<t>{n) dys 

^ r di/"^^ r 

+ V/ vCrUs-^da- (z/("Vw^)Cz/'^cia. (3.23) 

r^s=i Jd<t>{n) dyr Jd<t,{n) 

By passing to the hmit in fl3.23p as n — )• oo, by observing that XlfLi ~ ^ 
on d(l){Q) for all r = 1, . . . , iV, and by recalling that H = divz/ we obtain that 

/ VvVCu^dy + / vdiY{VCu^)dy = / vdivCda 

J(t>{n) J(P{n) JdcpiQ) 

+ [ Vv(^da- [ Hvu^da- [ ^(u'^da. (3.24) 

Jd^{n) Jd^(n) Jd<p{n) dv 



By combining fl3.20p and fl3.24p . we conclude. □ 

Proof of TheoremlQ We denote by R) the Banach space 

of the bounded real- valued bilinear symmetric forms in W^''^{Q) endowed with its 
usual norm. We denote by QiiW^'\n)f; M) the open subset oi Bs{iW^'^{n)f; M) 
of those bilinear forms which define a scalar product in W^''^{Q) topologically 
equivalent to the standard scalar product of W^'^{Q). Finally we denote by 
K[W^'^{Q)] ly^'^(n)) the Banach space of the compact linear operators in W^''^{fl) 
endowed with its usual norm. We note that the map from Aq to M) 
xK{W^'^{n); W^'^{n)) which takes G to the couple (< ■ >^,T^) is real- 
analytic, since it is the composition of real-analytic functions (see also the begin- 
ning of the proof of Lemma [3.13p . Moreover, by Lemma [3.121 is self-adjoint 
with respect to the scalar product < •, • >^ and the eigenvalues of T^f, coin- 
cide with the reciprocals of Aj[0]. Clearly, the set ^^[-F] coincides with the set 
{0 G : /ijM 7^ Vj G F, / G N\F}. Thus, by applying dH Thm. 2.30], 

it follows that ^f7[-F] is an open set in C'^{Q ; M^) and the functions which take 
G An[F] to 

^fA<P]= E ^'n[4>]■■■^^M (3-25) 

jl,...JhGF 

ji<-<jh 

are real-analytic for all h = 1, . . . ,\F\. Since 

A r^l F,\F\~h[<j)] , . 

^F,hm = TTT' (3-26) 

^F,\F\[(P\ 

for all = 1, . . . , where r^o[0] = it follows that Ap^cj)] depends real- 
analytically on G ^^[-F]- 

It remains to prove formula (13. 4p . Let G 6q[F], Xf[(P] and {fzjzgF be as in 
the statement. We set ui = vi o (p for all / G F and we note that {ui}i(zf is an 
orthonormal basis in VF^'^(f2) for the eigenspace corresponding to the eigenvalue 
A^^[0] of the operator T^. By [11] Thm. 2.30], it follows that 

drp,,(0)[^] = a^-'^hC^'j^^) < <^TMH,^i >h (3.27) 

^ ^ l&F 

for all G .4.q[F], ip G (7^(0 ; M^). Thus we have to compute the summands in 
the right-hand side of fl3.27p . For the sake of clarity, we believe it is better to 
compute < dTt^f'?/'] [m/], Mm >0 for all l,m E F. 

To shorten our notation we set A = Ai?[0] and ( = '0o0(^^). Since T^pUi = X~^ui 
for all / G -F, we obtain 

< dT^[lp][Ui],Um >cp 

=< iI^-A^Y-^^odJ4^][Tiui]+diI^-A^)^-^\^]oJ4Trui],Um ></, 

= (/^ - A^) [(/^ - A^)(~i) o dJ4^][Tr n^] + d(/^ - A^)^-'\^] o J^Tim]] KJ 

= dJM[TTUi][um] - A-M(/^ - A^M[ui][u^] . (3.28) 

By standard calculus, by means of a change of variables and formula fl3.16p it 
follows that 
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dA4V']NM= / Vvi{VC + Ve)Vvldy- [ VviVv^dwCdy. (3.29) 

J(f){n) J<f>{n) 

Since Af ; = vi on and ^ = Ai;; on (90(^2) for all / G -F, by the Divergence 
Theorem it follows that 

J(j>{n) 

~^J^{n) dysdyr dys dyrdyr dy^ ^Jd^n) " ' ^Vr dys 



dvi dvm , ^ f ^ d'^vi dvra . ^ Ovi d'^Vm , ^ d\ dv, 



rill rill f ^ I , riii riii riii riii riii-^ 



m 



9yr dys i<^(Q) dyrdys dy^ dyr dyl dyf. dy 

dyr dyrdys Jq^^q) \ dv dv J 

+ / divCVWzVw^ - {/^VrryVi + AviVVm)Cdy = \ [ {v^VVi 

J(h(n) Jd(h(n) 



J 4>(n) J^{n) 

- I VviVvl{vC^)da = A / V{viVm)C^da + / {VviVv^ 
Jd<t>{n) Jd<t>(n) J<t>(n) 

+vivjdw(dy- I (yviVvl + viVm){i^C^)da. (3.30) 

d()>{n) 



Moreover, by fl3.16p and a change of variable it follows that 



d/0[?/'][fi][fm] = / viVmdwCdy- (3-31) 
J<t>{n) 



By fl3:29D - fl33TD we obtain 



d(/^ - A^) [^] [vi] [vm] = / [{VviVvl + viv^)v - \V[viv^)\Cda (3.32) 



'd<t>(n) 

which combined with fl3.28p and fl3.15p yields 



< dT^[tlj][ui],Um >cf> 



A-1 



A ( h HviVrr, ] - (Vf/Vf^ + VlUn 



vC^da .(3.33) 
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By equalities (E2SD, fl^:^ and flS:^ it follows that 



dAFA#] 

= {dTp^iF\-h[<P][i']^F,\F\[(t)] - TF,\F\-h[(p]dTF,\F\ 



h+l-2\F\, 



\F\ - 1 

F\-h-l 



A 



h+l-2\F\ 



h 



F\ 
F\ 



X 



■Af Yl < dTMH^^i ><t>= -A^-^M0] /) E < 

IF 



h- 



+ vt -\f[4>\ 



IdF 
,2\ 



+ Hvf 



Formula fl3.4p then follows by observing that 

d{vf 



\Vvi\^ + vf-XF 
for all / G F. 



+ Hvf = \VtVi\' + (1 - A^[0]i/ - XM)vf 



□ 



4 Overdetermined problems 

In this section we consider isovolumetric and isoperimetric domain perturbations. 
Namely, given a bounded domain Q in of class C^, we consider domain trans- 
formations G satisfying either the volume constraint 

Vol = const, (4.1) 

or the perimeter constraint 

Per = const, (4.2) 

where Vol 0(f2) denotes the N-dimensional Lebesgue measure of 0(fi) and Per 
the (N-l)-dimensional surface measure of d(j){fl). It is here natural to introduce 
the functionals V and P from Aq to M defined by 

= Vol0(fi), and P[0] = Per0(fi), (4.3) 

for all G Note that by changing variables in integrals 

V[(p]= [ |detV0|rfx, and P[(p] = [ |z/(V0)"^ | |detV0| c/a, (4.4) 
Jn Jan 

for all G We recall the following 

Definition 4.5 Let fl be a bounded domain in M.^ of class . Let be a real- 
valued differentiate map defined on an open subset of Aq. 

(i) We say that G Aq is a critical point for T with volume constraint Iji4-i\ ) 
^f 

kerdV^(0) C kerdJ^(0). (4.6) 
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(ii) We say that </> G Aq is a critical point for T with perimeter constraint 

kerdP(0) C kerdJ^(0). (4.7) 

As is well-known definitions (i) and (ii) are related to the problem of finding 
local extremal points for 

min J^[(f)] or max J-'[4>], 

V[ij>]=const l/[iji]=const 



and 



min J^[<f)] or max J-" [(f)], 

Pf(/>l=const Pf(/>l=const 



respectively. Indeed if is a local minimum or maximum of a function J-" under 
condition (14. ip (resp. (14. 2p ) then inclusion (14. 6 p (resp. (14. 7p ) holds (see, e.g., 
Deimling [5] § 26]). 

Statements (i) and (ii) in the following theorem provide a characterization of 
all critical points of functions Ap^h under volume or perimeter constraint. We 
note that the proof of statement (iii) in the following theorem makes use of the 
same argument of [12^ Prop. 2.30] which is included here for the convenience of 
the reader. 

Theorem 4.8 Let Q be a bounded domain in of class C"^ and F be a finite 
subset of N. Assume that G Oq[-F] and that the eigenvalues Xj[(f)] have the 
common value Xplff] for all j G F. Let {vi}i(zf be an orthornormal basis in 
W^''^{(f){Q)) of the eigenspace corresponding to Xf[4']- 

(i) The transformation cj) is a critical point for any of the functions Ap^h, h = 
1, . . . , with the volume constraint ( [^.ip if and only if there exists Ci G M 
such that 

J2\^TVi\' + i^->^F[m~>^U<P]>f = Ci, ondcPin). (4.9) 

leF 

(ii) The transformation (p is a critical point for any of the functions Ap^h, h = 
1, . . . , with the perimeter constraint ( [^.^[ ) if and only if there exists 
C2 G M such that 

J2\^TVi\^ + i^->^F[<P]H-\l[^])vf = C2H, ond^in). (4.10) 

leF 

(iii) If (f){Q) is a ball then conditions Ij^^W and ^J^JW are satisfied. 



Proof. First we prove (i) and (ii). As for the case of finite dimensional vector 
spaces, the fact that is a critical point for Ap^h under volume or perimeter con- 
straint respectively, is equivalent to the existence of the corresponding Lagrange 
multipliers which means that there exist ci, C2 G M such that dAp^h{(f)) = cidV{(j)) 
and dAp^h^cj)) = C2dP(0) respectively (see, e.g., Deimling [5, Thm. 26.1]). By 
(I3.15p . (I3.16p . (14. 4p . by changing variables in integrals and the Divergence Theo- 
rem it follows that 

dV{(l))[ip]= [ (^ o 0(-i))i/^da, and dP{(p)[ip] = [ if (^ o 0(-i))z/^rfa , 

(4.11) 
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for all ip G C^(n;M^). Thus by combining fl3.4p . fl4.1ip and by the arbitrary 
choice of we deduce the validity of (14. 9 p and (I4.10p for suitable constants 

In order to prove (iii) we proceed as in [12] . Without any loss of generality, we 
assume that 0(f2) is a ball centered at zero. By rotation invariance of the Laplace 
operator, {vi o A}i^p is an orthonormal basis of the eigenspace corresponding to 
Ai?[0] for all A G OAr(R). Here OAr(R) denotes the group of orthogonal linear 
transformations in M^. Since both and {vi o A}i(zf are orthonormal bases 

of the same space, it follows that 



1^1 
1=1 



oA 



\F\ 

E 

1=1 



'I 5 



(4.12) 



for all A G 0„(R). Thus Yl\=^i'^i radial. This, combined with the fact that 
in the case of a ball H is constant, implies the validity of (14.91) and (I4.1UI) for 
suitable constants Ci, C2. □ 



Theorem 14.81 suggests to consider the following overdetermined systems where 
A is an eigenvalue of multiplicity m and ui, . . . ,Um is an orthonormal basis in 
W^''^{fl) of the corresponding eigenspace 



(PI) 



Ul, 



I ET=i I'^TVil' + {1-XH- X')vf = const. 



in Q , 
in dil , 

on dfl , 



(4.13) 



(P2)<^ 



Aui = Ul, 
^ = 

I = 1, . . . ,m, 

[ ET=i\'^TVi\' + ii-\H-x^^ 



in Q , 
in dil , 

consti/, on dQ . 



(4.14) 



It would be interesting to characterize those bounded domains Q in of 
class such that systems (PI) or (P2) are satisfied and to know whether the 
existence of solutions to one of these systems implies that f2 is a ball. 
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